Pion Condensation in Baryonic Matter: from Sarma Phase to 
Larkin-Ovchinnikov-Fudde-Ferrell Phase 
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We investigated two pion condensed phases in the frame of the two flavor Nambu-Jona-Lasinio 
model at finite baryon density: the homogeneous and isotropic Sarma phase and inhomogeneous and 
anisotropic Larkin-Ovchinnikov-Fudde-Ferrell(LOFF) phase. At small isospin chemical potential ^j, 
the Sarma state is free from the Sarma instability and magnetic instability due to the strong coupling 
and large enough effective quark mass. At large fij, while the Sarma instability can be cured via 
fixing baryon density ub to be nonzero, its magnetic instability implies that the LOFF state is more 
favored than the Sarma state. In the intermediate fij region, the stable ground state is the Sarma 
state at higher ub and LOFF state at lower hb- 

PACS numbers: 11.30.Qc, 12.39.-x, 21.65.+f 



I. INTRODUCTION 

Recently, the study on Quantum Chromodynamics 
(QCD) phase structure is extended to finite isospin 
density [J. The physical motivation to study QCD at 
finite isospin density and the corresponding pion super- 
fluidity is related to the investigation of compact stars, 
isospin asymmetric nuclear matter and heavy ion colli- 
sions at intermediate energies. In early studies on dense 
nuclear matter and compact stars, it has been suggested 
that charged pions are condensed at sufficiently high 
isospin density]! H S Q . 

While the perturbation theory of QCD can well de- 
scribe the properties of new QCD phases at extremely 
high temperature and density, the study on the phase 
structure at moderate temperature and density depends 
on lattice QCD calculation and effective models with 
QCD symmetries. The lattice simulation at finite isospin 
chemical potential 6] shows that there is a phase transi- 
tion from normal phase to pion superfluidity phase at a 
critical isospin chemical potential which is about the pion 
mass in the vacuum, /ij ~ ttItt. The QCD phase struc- 
ture at finite isospin density is also investigated in low 
energy effective models, such as the Nambu-Jona-Lasinio 
(NJL) model applied to quarks|l|l|i[ll[ll[l3 which 
is simple but enables us to see directly how the dynamic 
mechanism of isospin symmetry breaking operates. In 
the frame of this model, it is analytically proved[l3 UM 
that the critical isospin chemical potential for pion super- 
fluidity is exactly the pion mass in the vacuum, /ij = m-^ , 
and this relation is independent of the model parame- 
ters and regularization scheme. In addition, near the 
phase transition point, the chiral and pion condensates 
calculated in this model are in good agreement with the 
lattice simulation [6j. However, up to now, most of the 



'Email address: hely04@mails.tsinghua.edu.cn 
tEmail address: jin-m@mail.tsinghua.edu.cn 
■l-Email address: zhuangpf@mail.tsinghua.edu.cn 



studies in lattice and model calculations at quark level 
is done at zero baryon chemical potential, and the con- 
ditions of charge neutrality and beta equilibrium which 
are required in physical systems such as neutron stars 
can not be achieved. It is therefore necessary to study 
the effect of finite baryon chemical potential on the pion 
condensation in the frame of NJL model at quark level. 
Recently, this issue has been discussed in chiral limit [ig| . 

In a pion superfluid at zero baryon chemical poten- 
tial, the quark and antiquark of a condensed pair have 
the same isospin chemical potential and in turn the same 
Fermi surface (Strictly speaking, the term "Fermi sur- 
face" is only meaningful in weakly coupled fermi gas). 
When a nonzero baryon chemical potential is turned on, 
it can be regarded as a Fermi surface mismatch between 
the quark and antiquark. The Cooper pairing between 
different fermions with mismatched Fermi surfaces was 
discussed many years ago. In early investigation of super- 
conductivity in an external magnetic field, Sarma found 
a spatially uniform state[l3 where there exist gapless 
fermion excitations. However, compared with the fully 
gapped BCS state, the gapless state is energetically un- 
favored, which is called Sarma instability [lj|. Recently, 
such a spatially uniform ground state prompted new in- 
terest theoretically due to the studies on fermion pairing 
with large mass difference |l5j. two-component fermionic 
atom gas with density difference 1 161 Il7| . isospin asym- 
metric nuclear matter with neutron-proton pairing 
and neutral color superconducting quark matter[ 
I21I |22 |. It is now generally accepted that the Sarma 
instability can be cured in some cases, such as the sys- 
tem with l ong -range interaction where charge neutrality 
is reauired [20ll2l| and the system with proper finite range 
interaction between the two species of fermions with large 
mass difference 23]. In the Sarma state, the dispersion re- 
lation of one branch of the quasi-particles has two zero 
points at momenta pi and p2 which behave as two ef- 
fective Fermi surfaces, and at these two points it needs 
no energy for quasi-particle excitations. Since the sys- 
tem contains both a superfluid Fermi liquid component 
in the region p < pi,p > P2 and a normal Fermi liquid 




component in the region pi < p < p2, the Sarma state 
is now also called breached pairing state or interior gap 
state. 

Unfortunately, it is found that the Sarma state suffers 
negative superfluid density ,241 ^^ negative Meissner mass 
squared|23, |2a| which indicates a magnetic instability, 
and the LOFF state[23 123 with spatially non-uniform 
condensate is more energetically favored than the Sarma 
state|23,|30,|31|- In the study of isospin asymmetric color 
superconducting quark matter, isospin asymmetric nu- 
clear matter and atomic fermion gas with density differ- 
ence, the LOFF phase is proposed to be t he g round state 
and widely investigated |3^ 133 . yj, |33, ISg- However, 
it is recently argued that, the Sarma phase will be the 
stable ground state of asymmetric fermion superfluid in 
the Bose-Einstein condensation(BEC) rcgion^l6...17..37j. 
In relativistic fermion superfluids, the same conclusion 
is obtained when the coupling is strong enough and the 
fermions are heavy enough _38J. While the Sarma state 
may not be realized in color superconductivity, since the 
effective quark mass is very small at large baryon density, 
it is probably stable in a pion condensate matter with 
nonzero baryon density, since the effective quark mass 
is large enough near the phase transition point of pion 
condensation, i.e., around the isospin chemical potential 
/i/ = ttItt. In this paper, we study the Sarma state and 
LOFF state in pion condensed matter at finite baryon 
density in the frame of NJL model. 



strong coupling and large quark mass, it is therefore the 
stable ground state. At large /^/, the Sarma phase suffers 
the Sarma instability and magnetic instability. While the 
Sarma instability can be cured via fixed nonzero baryon 
density, the magnetic instability implies that the LOFF 
phase is more favored than the Sarma phase at large ^j. 
The schematic phase diagram of pion condensation at 
nonzero baryon density we obtained here is similar to 
the phase diagram of cold polarized fermi gas proposed 
by Son and StephanovflT]. and the isospin chemical po- 
tential here plays the same role of the coupling in atomic 
fermi gas there. 

The paper is organized as follows. In Sectional we 
review the formalism of pion condensation in the NJL 
model at finite isospin chemical potential and baryon 
chemical potential. In section IIIII we investigate the 
property of pion condensation at zero baryon chemical 
potential and show that there is a BEC-BCS crossover 
when the isospin chemical potential increases. In section 
IIVI we evaluate the Sarma pion condensed phase and dis- 
cuss the gapless quasiparticle spectrum. In section^ we 
study the thermodynamic stability of the Sarma phase. 
In section I VII we calculate the response of the Sarma 
phase to an external electromagnetic current and isospin 
current and discuss the dynamical stability of the Sarma 
phase. In section IVIII we construct the LOFF pion con- 
densed phase and discuss when it is more favored than 
the Sarma phase. We summarize in section IVIIII 
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FIG. 1: The schematic phase diagram of pion condensation 
in the /i/ — ub plane. 



II. NJL MODEL AT FINITE in AND /is 

In this section, we review the formalism of NJL model 
at finite temperature and isospin and baryon chemical 
potentials |llj. Since the isospin chemical potential which 
triggers pion condensation is large, /i/ > tti^, we ne- 
glect the possibility of diquark condensation which is 
favored at large baryon chemical potential and small 
isospin chemical potential. The key quantity to describe 
the thermodynamics of a system at finite baryon and 
isospin densities is the partition function defined by 



Z{T,iii,iXB.V)=TYe-'^^"-^^'^-^''''\ 



(1) 



The qualitative phase diagram of pion condensation 
at finite baryon density we obtained in this paper is il- 
lustrated in Fig^ in the isospin chemical potential and 
baryon density plane. The pion condensation can exist 
when the baryon density is not very high, otherwise the 
system will be in normal phase without pion condensa- 
tion because of the too strong mismatch. The critical 
baryon density for the transition from pion condensation 
to normal phase increases with the isospin chemical po- 
tential. The pion condensation starts at the critical value 
jjfj = mj^ and is separated into two phases at another crit- 
ical value /ij. At small /i/, the Sarma phase is free from 
the Sarma instability and magnetic instability due to the 



where V is the volume of the system, (3 the inverse tem- 
perature /9 = 1/T, iiB and fij are the baryon and isospin 
chemical potentials, and B and I3 the conserved baryon 
number and isospin number operators. If we take only 
quark field ■0 as elementary field of the system, the op- 
erators can be expressed as 



B= - I d^x^p-/o^ , 



h^l I d'^ij^or'^P, (2) 



the factors 1/3 and 1/2 reflect the fact that 3 quarks make 
a baryon and quark's isospin quantum number is 1/2. In 
the imaginary time formalism of finite temperature field 



theory, the partition function can be represented as 
Z{T,fii,HB,V) = f [dij][d'iP][dA]e- fo d- 1 d'^{c+^i^-fo4-) ^ 

(3) 
where C is the Lagrangian density describing the sys- 
tem, and /(, the quark chemical potential matrix in flavor 
space /t = diagd^uj l^d) with the u and d quark chemical 
potentials 

M-B , M/ Mb M/ /.x 

M« = - + y, M.= --y. (4) 

For convenience, we also use the notation fi = ij-b/S in 
the following. 

The flavor SU{2) NJL Lagrangian density is defined 
as 

C^i> (j^A'a^ - mo) V + G [(^V-)' + (Vi^srV')'] (5) 

with scalar and pscudoscalar interactions corresponding 
to a and tt excitations. The Lagrangian density has the 
symmetry Ub{1)^ SUi{2)<^ SUaC^) corresponding to 
baryon number symmetry, isospin symmetry and chi- 
ral symmetry, respectively. However, at nonzero isospin 
chemical potential, the isospin symmetry SUi{2) breaks 
down to Ui{l) global symmetry with the generator I^ 
which is related to the condensation of charged pions. 
At zero baryon chemical potential, the Fermi surfaces of 
u{d) and aTiti-d{u) quarks coincide and hence the con- 
densate of u and anti-d quarks is favored at sufficiently 
high fii > and the condensate of d and anti-u quarks 
is favored at sufficiently high fij < 0. We introduce the 
chiral condensate. 



and the pion condensates. 



(6) 



v2 



{^pij^T-ip) — v2{di'-f5u) = TT 



T2 



(7) 



with T± = (ti ± iT2) I \f2. A nonzero condensate cr 
means spontaneous chiral symmetry breaking, and a 
nonzero condensate tt means spontaneous isospin sym- 
metry breaking. The phase factor B related to the con- 
densates 7r+ and 7r_ indicates the direction of the C//(l) 
symmetry breaking. If the ground state is uniform, 9 is 
a constant and it does not alter the physical result. For 
convenience we choose 6* = in the following. 

In mean field approximation the partition function is 
simplified as 



Z{T, 



>/,A^B, V) = J[drP][di,]e--lo^^-l '*-^^™/ (8) 



■ J„^ dr J d^'^C^f 



with the mean field Lagrangian density 

£ = ^ [t^^d^ -m + Mn- iAriTs] - G{a^ + 'K^)i!, (9) 



where m — ttiq — 2G(t is the effective quark mass and 
A — —2Gtt is the effective energy gap. From the mean 
field Lagrangian density, we can derive the inverse quark 
propagator matrix in flavor space as a function of quark 
momentum, 

^_i. X ^ / 7^p^ + fiulo - m -«75A 

V -*75A 7^p^ + fidla - m 

(10) 

Using the formula listed in Appendix A, the thermody- 
namic potential in mean field approximation, 

n^-—\nZ^G{a^+n^)-—\ndetS-^ (II) 

can be evaluated as a summation of four quasiparticle 
contributions, 



n = G{a' + ^') - 6E / |^ff(^^(p))' ( 



12) 



where cui are the dispersions of the quasiparticles, 

uji{p) = E^ + ^1, uj2{p) = E^- 11, 

uj3{p) = E+ + ^i, Ui{p)=E+-ii (13) 

with the definitions 



E 



A 



(Eff + A2 , 



Ef =Ep± fii/2, Ep - ^p^+m^, (14) 

and the function g{x) is defined as g{x) = x/2 -I- Tln(I -|- 
g-s/T-j^ The gap equations to determine the effective 
quark mass m and pion condensate A can be obtained 
by the minimum of the thermodynamic potential. 



dm ' (9A ' drn^ ' dA-^ 
From the first order derivatives, we have 



(15) 



m — Too 



I2G 



d^ 



p m 



(2^)3 Ep 



^ (I - /(u;i) - f{co2)) 



+ T^(l-/(^3)-,/(^4)) 



A = 



I2GA 
I 



d3p 



E 



(27r)3 
(I - /(a>3) - f{uJi)) 



-^ (I - /(^i) - /(C.2)) 
^A 



A 



(16) 



with the Fermi-Dirac distribution function f{x) = 
1/ (e^'^ -|- l). This group of gap equations is invariant 
under the transformation /i/ -^ —fJ'i, and we can only 
concentrate on the case ni > 0. Once fi is known, the 
thermodynamic functions such as the pressure p, the en- 
tropy density s, the charge number densities ub and n/, 
the flavor number densities Uu and n^, the energy density 



e, and the specific heat c can be obtained by thermody- 
namic relations. 

To investigate the quark propagation in the pion con- 
densed matter, we should derive the explicit form of the 
quark propagator at mean field level 



exp li 



^ ' \ Sdu{p) Sdd{p) J 



(17) 



with the elements 



^7/,?A 



■>dd 



{po + fi + Ep) A+70 (po + M - E+) A_7o 



Sud 



^du 



{PO - UJ2){po +Wl) 

{pq + fJ- - E-) A^jQ 

{po - UJ2){P0 + ^l) 

-?:aa+75 ^ 

(po - L^2){pa + tJl) 
-iAA_75 



{po - UJ4){po + UJ3) 

{po + /i + E+) A+70 

{pO -UJ4){po +UJ3) 

-iAA_75 



(po -UJ4)iP0 +W3)' 

-iAA+75 



(po ~ uj2)iPo + LUi) ipo-LU4){pa+u;3)' 
where A± are the energy projectors 



A±(p) = i(l± 



70 (7 • p - 



En 



(18) 



(19) 



The occupation number of each kind of quarks, namely 
u, d, anti-u and anti-d quarks, is useful for our discus- 
sion in the following. They can be calculated from the 
positive and negative energy components of the diagonal 
propagators 5„„ and 5dd|l3i 



n+{p) - ulf{uj2)+vlf{-uJi), 

"d(P) = u'if{uJl)+V^_f{-UJ2), 

^diP) = u\f{LJ4)+vlf{~UJ3), 

IT-ZiP) = u\f{LJ3)+vlf{-UJ4), 



(20) 



in the notations n^ ^ stand 



where the symbols -I- and 

for quarks and antiquarks, respectively, and the coherent 

coefficients Wj_ and Uj. are defined as 



-l^^ 



E^ 



E^, 



^1m 



%y <- 



The number density of each kind of quarks and the total 
baryon number density are 



"uA 



6 



UB = 2 






^td(P)' 



(22) 



[/(u;2)-/(c.i) + /(c.4)-/(c.3)]. 



and zero baryon density, the QCD matter is a pure me- 
son matter, i.e., a Bose- Einstein condensate of charged 
pions. At ultra high isospin density, the matter turns to 
be a Fermi liquid with quark-antiquark cooper pairing[l|. 
Therefore, there should be a BEC to BCS crossover when 
the isospin chemical potential increases. In this section 
we discuss the gap equations (|16|) at T = ^.b = which 
has been investigated in 11] and point out that there are 
some signals of BEC to BCS crossover in the NJL model 
at finite isospin density. 

Since the NJL model is non-renormalizable, we should 
employ a hard three momentum cutoff A to regularize the 
gap equations. In the following numerical calculations, 
we take the current quark mass mo — 5 MeV, the cou- 
pling constant G = 4.93 GeV~^ and the cutof A — 653 
MeV. This group of parameters ensures the pion mass 
TOtt = 138 MeV and the pion decay constant f-^ — 93 
MeV in the vacuum. Numerical solution of the gap equa- 
tions is shown in Fig[21 For ^/ < ttt.^ the ground state 
is the same as the vacuum and the isospin density keeps 
zero, while for /z/ > m^ the pion condensate and isospin 
density become nonzero. The critical isospin chemical 
potential /Zj is exactly the vacuum pion mass m,r, and 
the phase transition is of second order. Since the chi- 
ral symmetry is spontaneously broken at small /i/ and 
almost restored at large /i/, the effective quark mass m 
plays an important role at small /i/ but can be neglected 
at large ///. In fact, the effective quark mass as a function 
of /!/ at /i/ > ttItt can be well described bvpT| 



to(0) ~ cr(0) 



m^ 
M/ 



(23) 




FIG. 2: The efFective quark mass m and the pion condensate 
A as functions of /ij/tti^. 



To explain why there will be a BEC to BCS crossover, 
we focus firstly on the dispersion relations of the 
fermionic excitations, 



III. BEC-BCS CROSSOVER AT ZERO /is 

It has been argued both in effective theory and lattice 
simulation that at finite but not very large isospin density 



(^1,2 



(p) = y(Vp2 + m2-/i,/2)' + A2, 

(P) = ]I{Vp^ + m^+ iJ^il'^ + A2. (24) 



In the case of /x/ > 0, W3 = L04 is the anti-particle excita- 
tion in the sense of isospin and irrelevant for our discus- 
sion. In Fig|31we showed the dispersion uji = uj2 at ^j — 
150, 450, 750 MeV. Obviously, the fcrmionic excitations 
are always gapped. At small /i/ with /i//2 < jn{fij), the 
minimum of the dispersion is at p = where the energy 
gap is \/m^ + A^ with /ijv = A'//2 — "m- which can be 
regarded as the corresponding non-relativistic chemical 
potential we will explain below. However, at large /// 
with fJ-i/2 > m(fj,j), the minimum of the dispersion is 
shifted to p ~ /^//S where the energy gap is A. This 
phenomenon is quite similar to the study of nucleon dis- 
persion at finite isospin density in the nonlinear sigma 
model|33. Such a phenomenon is a signal of BEC-BCS 
crossover. 
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FIG. 3: The dispersion relation of the fermionic excitation 
wi = W2 at p,i = 150, 450, 750 MeV. 
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FIG. 4; The occupation number of the fermionic excitation 
n+ = n- at ^/ = 150, 450, 750 MeV. 



Secondly, we discuss the occupation number of each 
kind of quarks , 






(25) 



In Fig^ we show the occupation number n+ ~ n^ at 
/i/ == 150, 450, 750 MeV. At small /z/ such as /i/ = 150 
MeV the occupation number is very small and smooth in 



the whole momentum region, while at large ni the oc- 
cupation number becomes large near p = and drops 
down with increasing momentum. When /i/ becomes 
very large, the occupation number is indeed of BCS type. 
This is the second signal of BEC-BCS crossover. 

Finally, we should point out some similarity between 
such a BEC-BCS crossover and the non-relativistic sys- 
tem. For small /z/ where m is large, the minimum of the 
dispersion is located at p = 0. Near this minimum with 
p ^ m, we can do the non-relativistic expansion 



V^T' 



p 

2m 



(26) 



and the corresponding non-relativistic chemical potential 
is 



fiN = — -m. 



(27) 



It is well known that for a non-relativistic Fermi gas, 
the chemical potential will become negative in the BEC 
region. Therefore, we can estimate the BEC region in 
our model by requiring /ijv < 0, 



-^ < m{fii) = m(0) —^ 
2 \ni 



(28) 



which determines another critical isospin chemical poten- 
tial /ij shown in the phase diagram in the introduction. 



fi° = [2m(0)m2; 



1/3 



(29) 



With the above chosen parameter set, we have ^j = 230 
MeV. When the effective quark mass in vacuum varies 
from 300 MeV to 500 MeV, /ij changes in the region 
230-270 MeV. Such a non-relativistic BEC in relativistic 
fermion superfluid was investigated in detail in J43 . The 
critical value /Xj we obtained here is similar to the one 

fi^j ^ (m„m-^) obtained in nonlinear sigma model|39| 
where to„ is the nucleon mass. Even though the above 
discussion is only qualitative, it will be important for us 
to understand the phase structure when a nonzero baryon 
density is turned on. 

The critical value f/j has another physical meaning. 
The chiral perturbation theory predicts that the chiral 
and pion condensates can be related to each other by a 
chiral rotation[l|, i.e., for /i/ > m^ there are 

cr(^/) = o-(O) cos a, 7r(/i/) = o-(O) sina (30) 

with cos a — m'^/ij.'j. We found that this behavior of 
chiral rotation holds approximately in the region ttt.^ < 
/i/ < fJj and breaks down for fij > ^^ in the NJL model. 



IV. SARMA PHASE AT FINITE fis 

When a baryon chemical potential is turned on, we 
may expect that the pion condensate will jump to zero 



at a critical baryon chemical potential which indicates 
a first order phase transition, like the behavior of chiral 
condensate. Such an assumption has been used in the 
studies of QCD phase diagram in the NJL model|9lll3|. 
However, when a possible Sarma phase is taken into ac- 
count, this will not be true. To see why a Sarma phase 
can appear, let us compare the dispersion relations of 
the fermionic excitations Wi, i = 1,2,3,4 obtained here 
with those obtained in the study of color superconduc- 
tor with charge neutralitv[l9|. The main differences be- 
tween the two cases are 1) the averaged Fermi surface of 
the two paired fermions is controlled by baryon chemical 
potential in color superconductivity but by isospin chem- 
ical potential in pion superfluidity, and the mismatch is 
served by isospin chemical potential in color supercon- 
ductivity but by baryon chemical potential in pion su- 
perfluidity, and 2) the effective quark mass is very small 
and can be neglected in color superconductivity but plays 
an important role in pion superfluidity when the isospin 
chemical potential is not too large. Taking into account 
these two differences, we can obtain from the conclusion 
for gapless color superconductivity|23| that, the Sarma 
phase can be realized in pion condensed matter when the 
gap A is smaller than one third of the baryon chemical 
potential, A < fi. 

From the possible solution of the gapless excitations, 
'^2{p) = and LU4{p) = 0, we derive the two zero points 
at momenta pi and p2 satisfying pi < P2, 



with 



pi = mJXl - 1, p2= m\ \l - 1 



Ai = 



A. 



(31) 



tii/2-^J^- 


-A2 


m 


/.,/2+Vm'- 


~A2 



(32) 



Obviously, for A > ^ or A < /Lt but |Ai^2| < 1, there is 
no real solution of the gapless excitation, all branches of 
quasiparticles are gapped. In this case there are 






(33) 



for all p, and ub = 0. Only in the case with A < /i there 
is the possibility to realize the Sarma phase. According 
to the values of Ai_2, there are three possible types of 
Sarma state. 

Type 1: Ai > and |Ai^2| > 1- In this case, only 
the branch W2 has two gapless nodes pi and p2, the other 
branches are all gapped. This case is similar to the gap- 
less two flavor color superconductivity z(j|. In this phase 
we have 



nt{p) 
"d(p) 



"d(p) 
Kip) 



(34) 



for p < pi and p > p2 and 

"« (P) = 1: "db)=0, 

njip) ^K{p) ^vl (35) 

for pi < p < p2. The baryon number density is 



hb 



pI-pI (Ai- 1)3/2 -(A? -1)3/2 



37r2 



37r2 



(36) 



the nonzero baryon number comes from the normal com- 
ponent in the region pi < p < p2- 

Type 2: |Ai| < 1 and IA2I > 1. In this case, only the 
branch a;2 has one gapless node p2, the other branches 
are all gapped. The system is in the paired state (|34|l at 
p > P2 and the state l|35|l without pairing at < p < p2- 
The baryon number density becomes 



Ub 



pI (Ai - 1)3/2 



37r2 



37r2 



(37) 



Type 3: Ai < and |Ai^2| > 1- In this case, the 
branch 102 has a gapless node p2 and the branch W4 has 
a gapless node pi, the other branches are all gapped. In 
this phase we have 



"ub) ^n^{p) 



(38) 



for p > p2 and 

n+(p) = l, n,(p) = (39) 

for < p < p2, and 

n+{p)^n-{p)^vl (40) 

for p > pi and 

nj{p) = l, n-{p) = (41) 

for < p < pi. The baryon number density is 



riB 



p?+pi_(Ai-i)^/^ + (A?-i)3/2 



37r2 



37r2 



(42) 



Before making numerical calculations, we firstly esti- 
mate where the Sarma phase is of type 2 only. From the 
judgement H32|l . it is clear that the condition to have only 
type 2 is m{p.i,^B) > /^//2. Since the effective quark 
mass at finite baryon density can not be larger than its 
value at zero baryon density, we have 



fii < ^J-i, 



(43) 



where /i^ is the critical isospin chemical potential we ob- 
tained in last section. 

Numerically solving the coupled gap equations ifTC)) for 
m and A, we can find all possible solutions, including the 
gapped and gapless (Sarma) states. At a fixed isospin 
chemical potential /x/, there is always a gapped solu- 
tion A(/i/,/iB) = A{^i,fiB = 0) = Aq provided ^ is 
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FIG. 5: The effective quark mass m and baryon chemical 
potential [i as functions of the scaled pion condensate A/Ao 
at /ij = 150, 450, 750 MeV in the gapless region. 



FIG. 6: The dimensionless quantities Ai and A2 determin- 
ing the type of Sarma phase as functions of the scaled pion 
condensate A/Ao at \ii — 150, 450, 750 MeV in the gapless 
region. 



less than the energy gap. Different from the result ob- 
tained in p3| where the gapless state can survive only in 
a narrow isospin chemical potential region, we found that 
the Sarma state can exist at any /x/ > TOtt- Obviously, 
the gap A in Sarma state is smaller than the gap Aq in 
gapped state due to the Fermi surface mismatch. 

The effective quark mass m and gap A in Sarma state 
are shown in Fig[Sl at different isospin chemical poten- 
tials. To show only the Sarma state, we have restricted 
the condensate in the region < A/Aq < 1. For /x/ < /x° 
such as Hi — 150 MeV, the system is still in chiral break- 
ing phase with large effective quark mass, and the baryon 
chemical potential /i starts at a value larger than Aq 
and then decreases slowly with increasing A in the whole 
gapless region. The behavior that fj, decreases with A 
is important, since it is directly related to the stability 
of the Sarma state, like the case in the BEG region of 
asymmetric fermion superfluidsjlg. For /i/ > f/j such 



as /i/ — 450 and 750 MeV, the chiral symmetry is al- 
most restored with very small effective quark mass, and 
the baryon chemical potential fi starts at a value smaller 
than Aq and then goes up in the gapless region. 

The quantities Ai and A2 which determine the type of 
the Sarma state are shown in Fig|Sl as functions of the 
scaled condensate A/Aq at fixed isospin chemical poten- 
tial. The gapless solution is of type 2 at /z/ < ^/, as we 
analytically discussed above, and type 1 at large enough 
fij. For intermediate /17, there exist both type 1 at large 
A and type 2 at small A. We never found a solution of 
type 3. In fact, this type of Sarma phase appears only in 
chiral limit where the effective quark mass is always zero 
in pion condensation region |l3|. 

Because of the Fermi surface mismatch in the gap- 
less state, the baryon number density of the system is 
nonzero. In Fig[7| we showed the baryon density ns 
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as a function of A/Aq. It is easy to understand that 
ub drops down monotonously with increasing A, namely 
with decreasing degree of mismatch, and goes up with 
increasing isospin chemical potential, namely with in- 
creasing average Fermi surface. Note that the baryon 
density ns and the ratio of ns to the isospin density 
nj here correspond, respectively, to the density differ- 
ence and the density asymmetry a in asymmetric fermion 
superfluidslla, IsJ Ell, and the ratio will decrease with 
increasing isospin chemical potential. 




FIG. 7: The baryon number density ub in the gapless phase, 
scaled by no = 0.16/m~^, as a function of A/Aq at fii — 
150, 450, 7507^61^. 



potential as a function of the condensate, the minimum 
corresponds really to the Sarma state. As is well known, 
to have a stable Sarma state at fixed chemical potential, 
the system should be under some specific conditions, for 
instance, a proper finite range interaction between the 
two species of fermions with large mass difference. Here 
we found a stable Sarma state at fixed chemical potential 
without mass difference in the frame of a four-fermion 
point interaction. Such a phenomenon was observed in 
a strong coupling modeljSq and is consistent with the 
argument in,16- I7j. 




A/A„ 



FIG. 8: The scaled susceptibility ka/m/ as a function of the 
scaled pion condensate A/Ao at /i/ — 150, 450, 750 MeV. 



V. THERMODYNAMIC STABILITY 

To ensure that the Sarma phase we observed in the 
pion condensed matter is stable, we should check both 
the thermodynamic and dynamical stabilities of the 
systemjl^ [iy. In this section we discuss the thermo- 
dynamic stability of the Sarma phase, i.e., the stability 
against a small fluctuation of the order parameter A. 
To this end, we expand the thermodynamic potential in 
powers of a small fluctuation 6 



r2(A + 5) - r2(A) 






(44) 



The linear term vanishes automatically due to the gap 
equation dfl/dA = 0. At zero temperature, the quantity 
'^A = hws^ which reflects the thermodynamic stability 
can be directly evaluated as 



KA = 



(2^ 



A2 



{Ezr 



V E^ 



S{UJ2) + (m/ 



-fJ-i) 

(45) 

where 0{x) is the step function of x. In FiglHlwe showed 
the scaled susceptibility ka/ij-'J in the Sarma state as a 
function of the scaled pion condensate A/Aq at several 
values of /i/. At /i/ < /ij, ka is always positive in the 
whole gapless region, which means that the Sarma state 
is free from the Sarma instability. In Fig|3 we plotted 
the thermodynamic potential at flxed baryon chemical 




A (MeV) 

FIG. 9: The thermodynamic potential n as a function of the 
pion condensate A at fixed chemical potentials fi = 260 MeV 
and m = 150 MeV. 

For /i/ > fJ,% the susceptibility becomes negative and 
the Sarma state suffers from the thermodynamic insta- 
bility. The question is, can we cure the Sarma instability 
at large fij by considering baryon number conservation 
in physical systems? When the baryon number is fixed 
to be nonzero which is required by charge neutrality, the 
gapped phase is ruled out, and the possible homogeneous 
and isotropic ground states are only normal phase with- 
out pion condensation and Sarma phase. In this case, 
the Sarma state may become stable like the gapless two 
fiavor color superconductivity |2£| . Different from the sys- 
tem at fixed baryon chemical potential where the ground 
state is determined by the thermodynamic potential Q, 
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FIG. 10: The thermodynamic potential at fixed baryon chem- 
ical potential and the free energy at fixed baryon number den- 
sity as functions of the pion condensate. 



the free energy F = 17 -t- fiEns controls the ground state 
for the system at fixed baryon number density ub- As a 
comparison, we plotted in Fig^|the thermodynamic po- 
tential at a fixed baryon chemical potential ji = 220 MeV 
and the free energy F at the corresponding baryon num- 
ber density ub = 2.62no in the Sarma state. While the 
Sarma state corresponds to the maximum of the thermo- 
dynamic potential and is therefore unstable in the case at 
fixed chemical potential, the minimum of the free energy 
is located at the Sarma solution and the Sarma state is 
really stable in the case at fixed baryon number density. 



VI. DYNAMICAL STABILITY 



We discuss now the dynamical stability of the Sarma 
phase, i.e., the stability against the perturbation of an 
external current. Since a non-zero expectation value 
(uij^d) breaks both electromagnetic f/em(l) gauge sym- 
metry and the global C//(l) symmetry, we will discuss the 
responses of the system to the external electromagnetic 
current J^^ — e-tpQ^^^ip with Q — diag{Qu, Qd) and to 
the isospin current J^ = ^'tpr^j'^tp. The former is just 
the electromagnetic Meissner effect, and the later corre- 
sponds to the superfluid density which is important to 
judge the stability of the Sarma phase |24l|. 



A. Response to Electromagnetic Current 

The quantity which describes the response of the sys- 
tem to an external electromagnetic current is the photon 
self-energy or photon polarization tensor defined by 



^Zip) = -2 






rTr 



r':^s{k)t:^s{k~p) (46) 



with f ^^ = eQj'^ at one loop level. In the following 
we discuss only the long- wave and static property of the 
response. The quantities that can describe this property 
are the Debye mass mjj and Meissner mass ttim defined 
as 



limn:^J^(a; = 0,p) 



1 



'"A/ - - o l™n(5y + P'^P3)KLi^ = 0, p). (47) 

Since there is no possible instability associated with the 
Debye mass, we focus on the Meissner mass only. After 
some algebras, the Meissner mass can be expressed in 
terms of the matrix elements of the propagator 5, 



'■M 



~2 / ii \-i \^-uX^u'VQ^dd^''^QuQdXud\ 



(48) 



with the definition 



n 

Tdd = T y^ Tr [SddliSddlt] , 

n 

Tud = T y^ Tr [SudliSduli] , (49) 

n 

where we should take the summation over the index i — 
1,2,3. To avoid the UV divergence, we should subtract 
the vacuum contribution, which ensures that mj^.j is zero 
when there is no pion condensation. For mj^.j > 0, the 
system exhibits diamagnetic Meissner effect, otherwise 
the response is paramagnetic which means a magnetic 
instability 25] . 

Taking 7^„, Tdd and T^d calculated in Appendix IBI we 
showed in Fig^Jthe Meissner mass squared mj^j for the 
gapless phase at fii — 150, 450, 750 MeV. We found that 
mj^j is positive in the whole gapless region for /i/ < //j 
which means stable Sarma state, and becomes negative 
in the whole region for large enough fij such as /i/ = 750 
MeV which indicates dynamical instability of the Sarma 
state. At intermediate fii such as /x/ = 450 MeV, the 
Sarma state is stable for small gap or large baryon density 
and unstable for large gap or small baryon density. 



B. Response to Isospin Current 

The quantity which describes the response of the sys- 
tem to an external isospin current is the isospin current- 
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FIG. 11: The electromagnetic Meissner mass squared in the 
gapless state, scaled by e^/if , as a function of the scaled pion 
condensate A/Ao at ^i = 150, 450, 750 MeV. 



FIG. 12: The mass squared mj in the gapless phase, scaled 
by /if, as a function of the scaled pion condensate A/Ao at 
m = 150, 450, 750 MeV. 



current correlation tensor defined by 



nrb) - -^ / (04 Tr [nS{k)T^,S{k - p)] 



(50) 



with Fg — 



5X37'^. We still focus on the long- wave and 
static property of the response by considering the quan- 
tity mj defined as 

mj = -- lim (<7y + p^Pj)^i{^ = 0, p) (51) 

which is proportional to the superfluid density ps of the 
system. After some algebras, mf can be expressed in 
terms of the matrix elements of the propagator 5, 



(27r) 



P (. 

7r)3 V 



% 



dd 



2%,. 



(52) 



Similarly, we should also subtract the vacuum contribu- 
tion to ensure zero m?i in the normal phase. 

We calculated numerically the mass squared m?j for the 
gapless phase at /x/ = 150, 450, 750 MeV and showed 
the result in Fig^l The pi dependence of the stability 
analysis of the Sarma state is very similar to what shown 
in Fig^Jfor the response to the electromagnetic current. 

In FigsE] and El the negative m\[ and mj behav- 
ior as /i(/i^ — A^)^^/^ and become divergent in the 
limit of A/Ao ~> 1 where /z ^ A, as observed in 
gapless color superconductivity [2^ and breached pair- 
ing superfluidity [2g. This divergence does not appear 
at small /i/ < /ij, since in the limit of /i — > A, the 
gapless solution does not appear due to the fact that 
the zero points of the excitations pi^2 become imaginary 
when M7/2 < m. This behavior is quite similar to the 
model [331 proposed recently. 

In conclusion, the Sarma pion condensed phase is al- 
ways thermodynamically stable at fixed baryon number 
density, and dynamically stable at small isospin chemi- 
cal potential fij < p'j but dynamically unstable at large 
enough /i/. For intermediate pi above and close to p.'j, 



the Sarma state is stable at large baryon density and 
unstable at small baryon density. 



VII. LOFF PION CONDENSED PHASE 

The order parameters Tr"*" and tt" are complex conju- 
gate to each other and can be set to be real only in a 
homogeneous and isotropic superfluid. In this section we 
consider the possibility that the phase factor 6 in the or- 
der parameter is non-uniform. Similar to the single plane 
wave LOFF state, we take the following ansatz for the 
pion condensate 



(V'«75T+V') = V 2(m75(i) = 7r+ = — „ 
{iPij5T_4>) = V2{di-f5u) = TT- = ^e-2«i-^. (53) 

We call this phase with a nonzero wave vector q 
the LOFF pion condensed phase. Note that this 
phase possesses the same symmetry of the p-wave pion 
condensation[3, |3j 01 and is continued with the LOFF 
phase defined at ultra high fir^. This phase is also sim- 
ilar to the chiral crystal phase|43 . After a transformation 
of the quark fields for each flavors. 



X„(x) = u(x)e ^'J-^, Xd(x) = d(x)e 



iqx 



(54) 



the Lagrangian with the new fermion fields in mean field 
approximation reads 

C = X [il'^df, - 771 + /i7o - T37 • q - iAri75] x 

-G(a2+7r2), (55) 

from which the inverse quark propagator matrix in the 
fiavor space as a function of quark momentum can be 
derived directly. 



S-\p,q)^ (56) 

7^Pa' - 7 • q + /"u7o - rn -«75A 

-J75A j'^p^ + 7 • q + fidla - m 
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and the thermodynamic potential is then expressed as 



n = G{a^ + TT^) 



T 

V 



In det S 



(57) 



To evaluate the thermodynamic potential, we employ the 
following formula which can be derived from the formula 
listed in Appendix EI 



[det 5-1 



1/2 



(po+M + e-)'-(e+-M//2)'-A' 



2A^ 



{po + fi- £-)' - (e+ + ni/2f - A' 



- q^ - VIp + qp + to^^Ip - q|2 + m^ 



(58) 



where e± are defined as 



e± 



2 (Vlp + qp 



±^p^ 



(59) 



The gap equations to determine simultaneously the 
condensates m and A and the pair momentum q are de- 
rived from the minimum of the thermodynamic potential, 
namely the equations 115() and 



dq 



d^n 



>o. 



(60) 



Numerically solving the gap equations is not a easy task 
due to the last term in H58(l . In the study of color super- 
conductivity, people simply dropped the last term since 
the gap there is small but the baryon chemical potential 
is large M^. Here it seems not good to employ this ap- 
proximation unless fj,i is very large. In this paper we will 
not directly calculate the LOFF state, but explain why 
the LOFF pion condensed phase is favored at large ^/. 
Following the treatment in ^3j ^M ■• ^^ expand the ther- 
modynamic potential in powers of the pair momentum q 
in the vicinity of q = 0, 



n{q) - 17(0) 



dq 



q=0 



I92f] 
2'df 



9=0 



(61) 



The linear term vanishes automatically because of the 
gap equation (|60|l . Separating the pair momentum de- 
pendence of the inhomogeneous quark propagator from 
the homogenous one. 



S ^(p,q)=5 ^(p,q = 0)-r37-q, 



(62) 



and taking the derivative expansion, we can easily obtain 



'df' 



(63) 



Obviously, for mj > at small /i/, the LOFF phase 
is less favored than the Sarma phase, and for mj < 
at large /j,/, the LOFF phase is more favored than the 
Sarma phase. 

From the calculation on the pion superfluidity at fixed 
baryon density in Section llVl we can determine the phase 




FIG. 13: The phase diagram of pion superfluidity in the /ii — 
ub plane at zero temperature, calculated in the flavor SU(2) 
NJL model. 



transition line from the pion superfluidity phase to the 
normal phase in the /i/ — ns plane at flxed temperature, 
and then with the stability analysis discussed in Sections 
IVl and IVII and in this section, the line which separates 
the inhomogeneous and anisotropic LOFF state from the 
homogeneous and isotropic Sarma state inside the pion 
superfluidity can be further fixed. The phase diagram at 
zero temperature is shown in Fig ll3l which is very similar 
to the estimation shown in Fig^ 



VIII. SUMMARY 

As a first step to achieve a physical system with pion 
condensate, we should consider the effect of a nonzero 
baryon density. In this paper we have investigated the 
Sarma and LOFF pion condensed matter at finite baryon 
number density in the frame of Nambu-Jona-Lasinio 
model, and checked the thermodynamic and dynamical 
stability of the Sarma phase. 

The isospin chemical potential /x/ controls the averaged 
Fermi surface of the quark-antiquark condensate, and the 
baryon chemical potential hb offers a mismatch between 
the quark and antiquark. At a fixed Hi > ^I'j = Jn^, 
the magnitude of the pion condensate decreases with in- 
creasing baryon density, and finally vanishes at a criti- 
cal baryon density where the system undergoes a phase 
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transition from the pion condensation to the normal 
phase. The critical baryon density increases with increas- 
ing isospin chemical potential. Inside the pion condensed 
phase, the system goes from the Sarma state to the LOFF 
state when /i/ increases. At low /i/ < /ij where fjj ~ 250 
MeV is determined by the pion mass and quark mass in 
the vacuum, the coupling is strong and the effective quark 
mass is large enough, the Sarma state is therefore ther- 
modynamically and dynamically stable. At high enough 
/i/ ^ /i^, the effective quark mass becomes very small, 
while the Sarma instability can be avoided via fixing the 
baryon number density to be nonzero, the dynamical in- 
stability of the Sarma state is still there. In this case, the 
LOFF state is stable. In the intermediate region with /i/ 
being above and close to /i/, the stable ground state is 
LOFF state at lower baryon density and Sarma state 



at higher baryon density. The above summary of the 
phase structure of pion condensed matter is illustrated 
in FigUl 

Finally, we should point out that the diquark conden- 
sation in two color QCD at finite baryon density has the 
same behavior as the pion condensation in QCD with two 
or three colors at finite isospin density. This has been 
confirmed by both lattice simulation and low energy ef- 
fective theory studies [i^. We expect that the same story 
will happen in two color QCD at finite baryon and isospin 
chemical potentials. 
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APPENDIX A: EVALUATING detS'^ 

To evaluate the thermodynamic potential Q at mean field level, we should calculate the determinant of the inverse 
fermion propagator with off-diagonal elements 



^ ^^'^)-(, -Z75A 5,-1 (p,q) 



(Al) 



where 



Sa ^ (P, q) = bo + l-J-a)lO -J -Pa -ma 



(A2) 



is the inverse propagator for the quasiparticlcs a = 1, 2 with the masses mi, 7712, the chemical potentials /ii, ^2 and 
momenta Pi = P2 = P in isotropic states and Pi = p -I- q, P2 = p — q in LOFF state. Taking into account of the 
identity 



we have 



detS-' ^det{j''S-'j''), 



detS-\p,ci) ^ det'/' ( , ,, -?"P?"A^ iriPi +P2), - [m, +m,))^J,A 

^ ' V ('y (^1+^2)^- ("^i+"i2))«75A mi -pi -A'' 

= [{pI -ml + A'){pi -mi + A') - A'Hp, + p,)' - {m, + m2)2)]' 



(A3) 



(A4) 



with Pa = {pa + /ia, Pa)- In the simplest case with mi = m2 = ?Tt, pi = p2 = p and fii = — /i2 = /i, we recover the 
BCS type result. 



det5 



-i(p) = {pI - iV^T^ - fif - A^y (pI - (TF^T^ + ^i)' -A^y. 



(A5) 



APPENDIX B: EVALUATING Tun, Tm AND %,d 



We evaluate the quantities T^u, %id and T^d in this Appendix. Firstly we calculate the trace in the spin space. 
Using the relations 



Tr[A±(p)7VA±(p)7V]=2|j, 

Tr [A±{ph'YA±iph'Y] - 2 ( 3 - |1 ) , Tr [A±{ph'YA^iph'j'] = 2| 



£;2' 
p 



Tr [A±(p)7°7^A^(p)7°7^] = 2 (^3 - ^) , 






,2 



(Bl) 
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we obtain 



iM 



^uu'~ii^uu^i 



,2 (JC^„ + ^ - E+f 



^p2 {iu},,+^i + E-f 

'El [(zc.„ + ^^Y - {E-^YY ^ 'eI [{^u,, + ^.f - {E+Yf 
iujn + ^i + E~ iiUn + fi- E+ 

EiJ (zc.„ + ^,Y - (E^Y (*^n + m)' - (eXY ' 



4 3 



,p2 {tu;„ + ^, + E+Y 






4(3-^ 



iw„ + ^- Ep iuJn + ^1 + E^ 

El) [lUr, + m)2 - (i^A)' («^n + M)2 - (^i)' ' 



]^^<^" 



dliSdul'i 



-A2 



,p 



-A2 



£2 [(zc.„ + ^)2 - {E-^Y? El [(^c.„ + /i)2 - {Ei) 



212 



4(3-^ 



-iA 



-iA 



Secondly we calculate the fermion frequency summations. From the decomposition of the summation 

[iuJn + ^. + eY 



^E' 



[{^OJn+^lY~E^Y 
for any constants e and E, and the simple frequency summations 

f{E + /i) + f{E~y)~l 



= a+{E + eYb + 2ec 



n 

^-^ Ml, 

n 



1 

(iujn + y^Y - E^ ~ 

1 Id 

[{iun +iiY -E^Y ' ^'dE 
1 



2E 

f{E + ^) + f{E - /i) - 1 
2E 



f{E + ^,) + f{E-^i)-l 1 df{E + fi) 



{iujn + ^i + EY{i(^n + H- E) iE^ 2E dE 

with the definition f'{x) ~ df{x)/dx, we can express some of the summations in (|B2|) as 



^E 

n 

^E 

n 

^E 

n 

^E 

n 

^E 

n 

^E 



{iuj„ + fi + E-Y 



(iuJn + fi - EpY 



W-W- ^V_f [LOi) +U_f [UJ2), 



{iLUr^+fiY-iE^y? 

A2 



E-, 



W_-U_ — hM_/ (tJi)+W„/ (W2), 

-^a 

/(^l)+/(^2)-l , 2 2 rw 



iiiu„ + ^iY ~ (E^VY 



2 2 



£;; 



+ u'_«^[/'(a;i) + /'M] 



/(C.3) + /M-1 , 4W, 



(i^„+/,)2_(^+)2]2 

(iw„ + ^ - E+Y 



= WlUl— : ^ ■ hi'+/'(t^3) + ut/'('^4), 



El. 



{iiOn + f^Y - (ElY 

A2 



Consider the decomposition 

iuJn + H + ei iuJn + M — ^2 



= u%v%_ — r h u\f (wg) + w^/ (a;4), 

^A 

2 2 /('^a) + /('^4) - 1 , 2 2 I fit X , .// M 



^Er 



(itj„ + ^)2 - £;2 (itj„ + /i)2 - £;2 



^1 



d+{E^- e{)[E2 + £2)6 - {E^ - ei)/ - (^2 + (^2)9 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 
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for any constants ei, 62, Ei and E2, and the simple frequency summations 



d = 



I = 



9 = 



n 

^E 
^E 

■n. 

^E 



1 



1 



1 



iuJn + fl- El ioJn + ^JL- E2 Ei - E2 

1 1 1 



f{Ei + Ai) + f{Ei -fi)-l f{E2 +^l) + I{E2 - m) - 1" 



(iw„ + ^lY - El {icun + m)^ - £^| 



i^?-i?l 



2£;i 



2£;2 



1 



1 



f{E2 - m) fJEi -^i) 1 - /(gi + ^l) 

{iuJn + fi)^ ~ Ef iu;„ + fi - E2 El -El 2Ei{Ei-E2) 2Ei{Ei+E2)' 

1 1 _ fjEi - fi) f{E2 -^i) 1 - .f{E2 + ^i) 

luJn + fi - El {luj.^ + fi)^ - El El~El 2E2{E2-Ei) 2E2{E2 + Ei)' 



(B7) 



we can express the other summations in (jB2ll as 



^E 



iuJn + n + Ep icun + IJ.- E+ 

(icon + tif - {E-^f {lUJn + tif - {E+y 



f{u;2)-ficJ4) , {E^-Ep){E++E+) 



El -El 



{Eir-{E+r 






ficUl) + f{LU2) - 1 /(C^3) + /(C^4) - 1 



2E-, 



2E1, 



[EXY - {E-^Y 2E^iE^ - E+) 2E^{E^ + E+) 
/(^2) , /M , l-/(c^3) 



^{Ea ^ Ep 
'{E+ + E^ 



I{^2) 



l-f{^l) 



[E^Y-{E+Y 2E+{E+-E^) 2E+{E+ + E^) 
iuJn + fJ- — Ep iuj„ + fi + -B+ 



itiOr. + /X)2 - {E^Y {itOn + /i)2 - (E+Y 



Ea ^ E^ 



-{e^+e;) 

HE+-E+) 

^E 



iE^r (ElY 



fiuJi) + f{L02) - 1 .f{LO^) + /K) - 1 



2E-, 



2E1, 



.fi^2) 



1 - /(^l) 



[E+Y ~ {Ea? ^E^iE^ - Ei) 2E^iEl + E+) 
f{^2) , /M , l-/(c^3) 



(£;^)2 - {ElY 2E+{E+ - EZ) 2E+iE+ + E^) 
A A 



{lUn + m)' - {E-^? {l^n + fiY - (E+Y 

A2 \f{uJi)+f{uj2)-l f{uj3)+f{uJi)-l 



{e-^Y^{eX? 



2E 



A 



2Ei 



(B8) 



For a general combination 



T(a, (3) = oYTuu + ii'^Tdd ~ 2a(iTud, 



(B9) 



we can separate it into a diamagnetic part Td and a paramagnetic part 7^, 



r(a,/3) = 2^rp + 4(3-^)r. 



(BIO) 



with 



Tp = (a-pyutv 






2 \2 /•// 



+ {avi + PuiYfiiUi) + [aut + (3vtY fi^^) 



% 



+ {a~l3Yu\v 



[E-^Y [Elf 



1 2 2 /(^S) + /('^4) - 1 



El 



+ (m4 + /3«2 )2^,(^^) ^ (^2 ^ /3w2 )2/'(^4)^ 



2a/3A^ + (a^+/3^)((ii;^)^-i^-£;;) 



+(a^ - /?") 



2xEi^E+ - E^Ep 



(E^Y-iEXr 
a2-/32 



/(a;i) + /(c^2) - 1 
2E^ 



(/K) + /M - 1) - (m/ 



- (m/ 



(i?i)^-(i^zF 



i?p"/M + ^(i + |^)/M 



15 



(Bll) 



-M/J 



-A*/J 



^(l-|^Hl-/K)) + (/i/ 



-A^/) 



We set a — Qu = 2/3 and /3 — —Qd = 1/3 for calculating m|,j and a = /? = 1 for mj. Note that only for a ^ (3, the 
formula is invariant under the exchange fii — > — ///, otherwi se the isospin symmetry is explicitly broken. 
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